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We develop a method to compute shakeup effects on photoluminescence from a strong magnetic 
field induced two-dimensional Wigner crystal. Only localized holes are considered. Our method treats 
the lattice electrons and the tunneling electron on an equal footing, and uses a quantum-mechanical 
calculation of the collective modes that does not depend in any way on a harmonic approximation. We 
find that shakeup produces a series of sidebands that may be identified with maxima in the collective 
mode density of states, and definitively distinguishes the crystal state from a liquid state in the absence 
of electron-hole interaction. In the presence of electron-hole interaction, sidebands also appear in the 
liquid state coming from short-range density fluctuations around the hole. However, the sidebands in 
the liquid state and the crystal state have different qualitative behaviors. We also find a shift in the 
main luminescence peak, that is associated with lattice relaxation in the vicinity of a vacancy. The 
relationship of the shakeup spectrum with previous mean-field calculations is discussed. 



PACS numbers: 78.20.Ls, 72.20.Jv, 73.20.Dx 



I. INTRODUCTION 

It was first pointed out by Wigner^ more than sixty 
years ago that an electron gas will undergo a zero- 
temperature, quantum phase transition into a crystalline 
phase as the density is lowered, since quantum fluctua- 
tion effects diminish more rapidly than Coulomb correla- 
tion. Only two decades ago, the flrst convincing evidence 
of an electron crystal was presented for a system of elec- 
trons on a He surface^ . The electron densities attainable 
in this fashion are extremely low, however, making this 
an unattractive system for observing the quantum phase 
transition. Semiconductors are much more attractive sys- 
tems in this sense, because one has great control over the 
electron densities, through dopant concentrations. A par- 
ticularly good candidate for observing the Wigner crystal 
(WC) is the two-dimensional electron gas (2DEG), as re- 
alized in modulation doped semiconductors. Samples of 
this type are now available with such high quality that 
the electron groundstate is not necessarily dominated by 
disorder. The possibility of observing the WC is further 
enhanced by the application of a strong perpendicular 
magnetic fleld, which quenches the kinetic energy, and 
allows the formation of a crystal state at higher densities 
(for which disorder effects are less important) than would 
be possible without it. 

Experimental evidence which may be associated with 
the WC in 2DEG's has accumulated over the last sev- 
eral years^, including rf data, transport experiment, cy- 
clotron resonance, and photoluminescence (PL) experi- 
ments. In the last of these probes, which has produced 
much intriguing data, either a valence band hole'* or a 
hole bound to an acceptor^'® recombines with an elec- 
tron in the 2DEG, producing a characteristic photon 



spectrum. A mean-fleld analysis of the latter type of 
experiment showed that the PL spectrum has, in prin- 
ciple, characteristic signatures of the WC: a "Hofstadter 
butterfly"® spectrum for the case of weak interactions 
between the electrons and the hole, and a characteristic 
shift in the PL spectrum upon melting of the crystal. 

In this paper, we go beyond the mean-field approxi- 
mation, to examine shakeup effects on the PL spectrum; 
i.e., we will examine how the collective mode spectrum of 
the WC (which, at long wavelengths, corresponds to the 
classical phonon spectrum) , and the fact that some of 
these modes may be excited in the electron-hole recombi- 
nation process, modify the results of the mean-field the- 
ory. We will consider in detail only the case of a localized 
hole®'®. Our method is purely quantum-mechanical, and 
treats both the tunneling electron and the other lattice 
electrons on the same footing. Furthermore, we employ 
a quantum treatment of the collective modes to realisti- 
cally account for contributions both from small and large 
wavevector collective excitations of the lattice. This ap- 
proach has the further advantage of depending in no way 
upon a harmonic approximation for the lattice^. The 
method allows in principle for shakeup of arbitrary num- 
bers of these excitations. Since we are working in the 
strong magnetic field limit, we consider only excitations 
within the lowest Landau level (LLL). Our main results 
are: (1) Shakeup effects shift the main PL peak to higher 
energies than found in a mean-field treatment^. (2) The 
Hofstadter spectrum is eliminated from the PL spectrum, 
even in the case of weak electron-hole interactions® (al- 
though we will argue below that it survives in the itin- 
erant hole case^). The sudden shift of the PL spectrum 
upon melting, by contrast, survives even when shakeup 
is included. (3) Phonon sidebands appear that corre- 
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spond to maxima in the phonon density of states (DOS); 
some (but not all) of these sidebands are results of van 
Hove singularities in the DOS, and so are characteristic 
of an ordered WC state. For the case of weak electron- 
hole interactions®, we do not see these sidebands in the 
liquid state, so that phonon satellites uniquely distin- 
guish between a liquid and a solid state. Interestingly, 
for stronger electron-hole interactions, a shakeup satel- 
lite persists even above the melting temperature. We 
expect this sideband to lose oscillator strength relative 
to the main peak, either with increasing temperature or 
decreasing electron-hole interaction strength. The latter 
may be accomplished by examining PL from several sam- 
ples with different acceptor - 2DEG setback distances® . 

This paper is organized in the following way. In Sec 
II, we review the mean-field theory for the photolumi- 
nescence, as well as the time-dependent Hartree-Fock 
approximation for collective modes in strong magnetic 
fields, and show how they may be combined to calculate 
the phonon shake-up effect on photoluminescence. We 
present and discuss our numerical results in Sec. Ill, and 
make some concluding remarks in Sec. IV. A brief account 
of some of these results has been published previously. 

II. THEORY 

In this section, we first set up the theoretical frame- 
work by reviewing the general expressions for the time- 
dependent Hartree-Fock approximation (TDHFA) of the 
2DEG in a strong magnetic field and outlining the mean- 
field theory for photoluminescence from Wigner crystal. 
Then we introduce the expressions for collective modes 
in Wigner crystal. Finally we present the detailed for- 
malism to calculate the collective modes shake-up effects 
on the photoluminescence. 

A. Single-Particle Properties 

It is well known that non-interacting two-dimensional 
electrons in a perpendicular magnetic field (B = —Bz) 
have an energy spectrum of discrete Landau levels: iJjv = 
[N + \)nuJc,N = 0, 1,2, where = eB/m*c is the 
cyclotron resonance frequency. Working in the Landau 
gauge(A = —Bxy), and with periodic boundary condi- 
tions in the t/-direction, the single particle eigenstates are 
given by: 

<Y\NX>=—exp{iXy/ll)(t)N{x-X). (1) 

Ly 

Here Ic = (hc/eBY^^ is the magnetic length and (/^at is 
the one-dimensional harmonic oscillator eigenstate with 
oscillation centers X. The allowed values of A are sepa- 
rated by 2iTl'^lLy. The degeneracy of each Landau level 
is given hy g = S/^irl'^, with S the area of the 2DEG. 



Das Sarma Group Preprint, 1996 

To address the problem of photoluminescence from 
a strong magnetic field induced two-dimensional WC, 
we now consider a system with an interacting 2DEG 
and a layer with a low density of holes (either local- 
ized or itinerant) separated from the electron plane by 
a distance d, in the presence of a strong perpendicular 
magnetic field. We assume the magnetic field is very 
strong; t.e., Tiojc » e^/ea, and the electronic filling 
factor v = nhc/eB << 1 in the WC regime, so that 
we can assume only the lowest Landau level is occupied 
by electrons. The general many-body Hamiltonian can 
be expressed in terms of single-electron eigenstates (only 
N = eigenstates are included) as follows: 

H = Ho -\~ Tiee , (2) 

where 

X i 

'^ee = ^J2 J2 Vcid) <Xi\exp(i(i-r)\Xi> 

q7£0XiX2X3X4 

X < X2\exp{-iq ■ r)\X3 > a\^a\^ax-^ax^ 
'^-'^ = J2 J2 (q) < Ai \exp{i(i ■ r) | A2 > 

q X1X2 ij 

X < i\exp{-iq ■ r)\j > a^^ax^clcj. (3) 

Here is the single-particle zero-point energy which 
is constant. Wee is the electron-electron interaction, in 
which Vc{(\) = 2iTe^ jeq is the two-dimensional Fourier 
transform of the Coulomb interaction. "Heh is the 
electron-hole interaction, in which Veh{^) is the Fourier 
transform of the interaction between one electron and one 
hole. For a system with a small density of holes, it is a 
good approximation to also ignore hole-hole interactions, 
which are extremely small in comparison with "Heh- In 
our current consideration, we also do not include the ef- 
fect of impurities and any external potential. In practice, 
we may also drop because it is just a constant. The 
matrix elements in Eq.(3) are given by 

< Ai \exp(i(\ ■ r) | A2 > 

= exp(-q4Xi + A2) - ^)Sx,,x,+g,P^ (4) 

It is well known that the ground state configuration for 
a two-dimensional Wigner crystal is a triangular lattice^ ^ . 
In the presence of electron-hole interaction, the Wigner 
lattice should deform around the position of the hole. In 
order to take into account this effect, while still taking 
advantage of the periodicity of the WC, we divide the WC 
system into hexagon unit cells. Each supercell contains 
a finite number of electrons, and one hole (in the center, 
for localized hole case). The finite size of our unit cells 
will not be significant for large enough supercells. 
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Denoting h as the distance between the centers of nearest-neighbor supercells, then the basis vectors for the super- 
lattice are Li = &(1,0,0) and L2 = &(^,i,0). Then the basis vectors for the reciprocal lattice are given by 

Gi = Go(l,0,0) and G2 = Go(-^, ^,0), where Go = {j^bj \ 

In a lattice system, all the expectation values for single particle properties in momentum space are non-zero only 
for momenta equal to reciprocal lattice vectors. We first define the density operator: 

n[G) = / d'^rexp[—iG ■ r)n[r) = gp[G)exp[ 4^)' ('^) 

where 

P(^) = ^ E e-^«-(^^+^=)<5x„x.-G„/=at,^ax,. (6) 

Xi X2 

One may easily show that < p{G = 0) >=< Ng > / g = v, where Ng is the electron number operator. 
In terms of many-body operator decomposition, Hartree-Fock approximation can be expressed as: 

«k"k"-^3"-^4 =< a\^ax^ > a^,ax3- < dx^^Xs > a^x^^x, (7) 

In this effective mean-field theory, a single electron moves as if in an average external potential provide by other 
electrons through both the direct (first term) and exchange (second term) interactions. For localized holes which 
are very far away from each other, we can also apply the approximation that < cjcj >= Sij . After some algebra, 
Hartree-Fock Hamiltonian for the 2DEG (with an electron-hole interaction appropriate for a localized hole) in the 
lowest Landau level can be written as: 

=gJ2 [^(^) < P(^) > +nhVeh{G)e-^"'"/^] p{G) (8) 

G 

where rih is the density of holes. Note that we ignore the negative sign in front of G because of apparent symmetry 
in this problem. W{G) is the effective Hartree-Fock interaction: 



-G 



2]2 



(9) 



where Io{x) is the modified Bessel function of the first kind^^. 
We define single electron Green's function: 

G{Xi,X2;t) = -< TraxdT)a{^{0) > (10) 
It is convenient to define the Fourier transform 

G(G,r) = - V e-^«-(^^+^=),5x„x,+G„;=G(Xi,X2;r). (11) 

g Z / 

X1X2 

We will use this form of the Fourier transform throughout this work. The TDHFA is derived by writing the equation 
of motion for the Green's function. 



^G{X,,Xr,T) = -^ 

OT OT 

= -Sx.xAr)- < Tr[H - iJeNe,axAir)a^xM > 



G(Xi,X2;r) = -— < TraxAr)a^xM > (12) 



One can compute the commutators explicitly, and simplify the result using a Hartree-Fock decomposition presented 
in Eq.(7). After Fourier transforming with respect to time, the equation of motion for G(G,w„) can be written as^^ 

(ic^„ +//e)G(G,ic^„) -^5(G,G')G(G',ic^„) = Sg,o, (13) 

G' 

where 
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5(Gi, G2) = [W(Gi - G2) < p(Gi - G2) > 

+nnVen{G, - G2)e-(°-°=)='^H e'(°^ x°=)'^/^ (14) 



We can directly diagonalize matrix B and obtain its eigenvectors Vj{G) and eigenvalues w|, after which the Green's 
function can be written as 

Vj{G)V*{G = {]) 

^ lUJn + //e - 

The density of states for electrons is then given by: 



D{E) = --/m(G(G = 0,E+i<5))(|) 



Finally, the density operator can be expressed as 



<p{G) > =G(G,r = 0-) 

= J2^j(G)V*(G = 0)fpD(i^]-^le). (17) 

j 

Here fpoix) = [1 + exp{l3x)]~^ is the Fermi-Dirac distribution. Since < p{G = 0) >= v, we can self-consistently 
calculate the chemical potential pe, the density of states and the electron density in G space. By iteratively solving 
Eqs.(13),(15), and (17), we can calculate the density configurations for the Wigner crystal. 



B. Mean-Field Theory for Localized Holes 

We now present our theory for photoluminesence from the WC in a strong magnetic field. The photoluminescence 
intensity is given, for a single localized hole state, by 

P{^) = ^J2J2 e-^"/'"^ I < m, 0\L\n, h > \''S{l, - E„ + E^), (18) 

n m 

where Z = S^e""^"/''^"^, \n, h > is a many-body electron state with energy £"„ and N electrons when there is a core 
hole present, |m,0 > is a many-body electron state with N — I electrons and energy Em, 1^ is the luminescence 
frequency, and L = J d'^rip[r)iph{Tc) is the luminescence operator, with ij{x) the electron annihilation operator and 
tph{x) the hole annihilation operator. As written, the initial state is actually higher in energy than the final state, 
and we find it convenient to rework the problem in terms of absorption rather than emission. To accomplish this, we 
add a term H' = —EqCqCq to the Hamiltonian, where cj creates a localized hole, and take the limit Eq 00. It is 
not difficult to show 

P(^)= Km Elit^LllA (19) 

^ ' Bo^oo no(^o) 

where P' is the absorption spectrum of the new Hamiltonian, and no is the average occupation of the hole state, 
which just becomes one in the limit Eg — ?■ 00. The absorption spectrum is identical to Eq.(18), except one needs to 
add the energy Eq to all the quantities E„ in the expression. After standard manipulations^^, one can show that 

P'(.) = ^^^-^ImR(. + iS) (20) 

The function R{ui + iS) is a response function, which continued to imaginary frequency has the form 
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Tl{tbJn) 



(21) 



To compute this quantity, we consider (for the case of a localized hole state) instead of a single hole, a periodic 
(hexagonal) lattice of them, with a unit cell that contains as many electrons as can be handled numerically. We 
allow neither interactions between the holes nor tunneling between the hole sites, so that in the limit of large unit 
(super)cells, one should expect the result to be the same as for the isolated hole case. Expanding the electron and 
hole creation and annihilation operators in L in Eq.(21) in terms of electronic states in the lowest Landau level and 
localized hole states, one finds 



RH = S|E^(G,-)e-^^'^/^ 

where i?(G,w) is the Fourier transform of the following quantity: 

nij{Xi,X2;T) = -< TraxAr)ci{T)c]{0)a^x^{0) >, 



and 



X1X2 



We write down the equation of motion for lZij{Xi, X2', r) in terms of its commutator with the Hamiltonian: 



^TZ.^{X,,Xr,T) = -^< T,ax,(r)c,(r)c](0)at^^(0) > 
- < [axiCi,c]a^J > S{t)- < Tr[Heff - /u(A^e - #ft), axiCi](r)c]a^^ >, 



(22) 

(23) 
(24) 

(25) 



where Tieff = 7i — Eo'^-Cicj and Ng and are corresponding electron and hole number operator. Following 
steps closely analogous to those used to find the equation of motion for the Green's function in Sec. A above, the 
Hartree-Fock approximation for Eq.(25) becomes^'^^ 



(w + iS)Rij(G,Lj) = < p(G) > Sij + (Eo - -hcjc - Eh)Rij(G,Lj) 



(26) 



J2 WiG') < p{G') > e-'(°x°')''/2i?ij.(G - G', 

G' 

l^J2v,,i-G') < p{G') > e-«'='^/X-(G,c. 



2ttI 



nh 



J2 V;4G')e-«"'^/4-'(GxG')'^/2i?,^.(G - G', 



(27) 



It is apparent that the solution to this satisfies Rij{G,ui) = R{G,ui)Sij. This result may be expressed in the form 



E 

G' 



{iu + iS + Eo- iUo)SG.G' - B{G, G') 



R{G',ui) 



<piG) >, 



where 



Wo 



(28) 



(29) 



and B is exactly given by Eq.(14). 

We see that the form of R is essentially that of electron Green's function as in Eq.(13). By inverting Eq.(28), we 
have 



R{G,oj) 



■E 

G' 

■E 

jG' 

■E 



(w + iS + Eo- wo)<5g,g' - 5(G, G') 

V,{G)V-'{G')<p{G')> 
uj + iS — ujo — i^j 

We{G,j)fpD{i0l -p,) 

UJ + iS — ujo — w| 



< P(G') > 



(30) 
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Here we have already dropped Eg because it cancels out when we calculate the final photoluminescence power using 
Eqs.(19) and (20). We can see that R has poles at precisely the same energies as the poles in the electron Green's 
function for the system in the presence of the external interaction Veh due to the hole^'*, up to the constant energy shift 
ujo- However, the residues of the poles are not the same as for the Green's function. The residues are determined by 
the overlaps of the mean-field single particle electron wavefunctions with the hole wavefunction. It is thus appropriate 
to think of i? as a weighted Green's function. Since -P'(w) involves the imaginary part of R(w), the photoluminescence 
spectrum at this mean-field level represents a weighted measure of the single particle density of states of the electron 
system with a localized hole. 
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FIG. 1. Electron density of states and PL spectrum for perfect WC with no electron-hole interaction for = 2/7. (a) Electron 
DOS below chemical potential at zero temperature (inset: DOS above chemical potential), where Ee = \iOc', (b) PL spectrum 



at different temperatures: 
dash-dotted). 



T = 0.0045T„eit (solid), T = 0.45T„eit(dotted), T = 0.9T„eit (dashed) and T = 1.12T„eit(inset: 



The results of this mean-filed approach have discussed previously^, so we only briefiy summarize them here. The 
DOS for a perfect WC is a Hofstadter butterfiy. For fractional filling factor v = p/q, this has q subbands, and at 
zero temperature, the lowest p bands are filled. We should expect to observe this structure; i.e., p lines for v = p/q, 
in the photoluminescence spectrum of an ideal Wigner crystal if we turn off the electron-hole interaction. Presented 
in Fig.l are the density of states (a) and the photoluminescence spectrum (b) for an undisturbed Wigner crystal at 
v = 2/1 with no e-h interaction. We can clearly see that the DOS has 7 bands, with only the lowest 2 occupied at 
low temperature. As expected, the photoluminescence has 2 peaks with a splitting identical to that of the DOS. An 
observation of this behavior in experiments would directly confirm the presence of a Wigner crystal in the system. 
Furthermore, as a function of temperature, the PL spectrum shifts very suddenly as the crystal melts to higher 
frequencies. This energy shift represents the latent heat of melting of the crystal. Once again, observation of this 
behavior would directly confirm the presence of a WC in the system. 

While these effects are relatively simple to understand within the mean-field approach, one needs to question 
whether they can survive fiuctuation effects. In particular, the smallness of the Hofstadter gaps in Fig. 1 suggest 
that they may not be robust enough to survive beyond the mean-field level. We will see for localized holes that this 
is indeed the case, and that this structure is replaced by a spectrum that refiects the collective mode density of states 
rather than the mean-field single particle spectrum. To show this, we proceed by reviewing how the collective modes 
are computed in the TDHFA, and then show how they may be incorporated into the PL spectrum. 



C. Collective Modes in Wigner Crystal 



We define the dielectric function as the density response of a two-dimensional electron system to an external 
perturbation ?7(p,r), i.e. 



-x(p,(i;t-t 



(31) 



SU{p,T) 

In the special circumstance of Wigner crystal, the density-density correlation function has the following special form: 

x(Gi + q, G2 + q; T) = -g< Tp(Gi + q, r)p(-G2 - q; 0) >, (32) 
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where p = p— < p >. From now on, we denote x(Gi + q, G2 + q; r) as XGiGjIci; ''")• The Fourier transform of this 
response function contains poles at the collective mode frequencies of the system. This function can be numerically 
computed using a generalized random phase approximation (GRPA)^^ which does not assume small displacements 
of the lattice electrons (as is necessary in classical approaches^), and thus gives a realistic dispersion relation for 
the collective modes across the entire Brillouin zone. In this subsection, we present the formalism for this GRPA 

calculation including the electron-hole interaction. (The magnetic length Ic = [Tic/eBj will be set to unity). 
In the real space, the dielectric function is given by: 

x{Xi,X2;X3,Xi;T) = -g < Trp{Xi, X2){t)p{X3, Xi){0) >, (33) 
where p[Xi,X2) = a^^axj— < >• The equation of motion for x can be written as: 

dx{Xi,X2;X3,Xi;T) 
dr 

= -g <[p{X,,X2),p{X3,X4)]> S{t) 

-g <Tr[n- pN, p{Xi , X2)](r)p(X3, X4)(0) > 
= -Ti-T2, (34) 

where 

T, = g <[p{X,,X2),p{Xs,X^)]> 5{t) 

= 9[< dx^ax^ > Sx^x^ - < al^^ax^ > <^XiX4]<^(t)- (35) 
We define the Fourier transform of this quantity as 

7^i(q,q';r) = ^ ^ e-^^^(^^+^=),5x„x.-,„ (36) 
- J2 e-^^^(^^+^^)<5x3,x,+giTi(Xi,X2;X3,X4 



^ „ u;r 

g ^ — ' 

X3X4 



so that in frequency and momentum space, 

f'i(G + q,G' + q;w) =< p{G - G') > Jg-'f^+i) >< (G'+q)/^ 

_gi(G+q)x(G'+q)/2l (37>) 



The most important term T2 comes from the electron-electron interaction and electron-hole interaction in the Hamil- 
tonian. Following the method of Ref.l2, a Hartree-Fock decomposition is applied to this term, to obtain 

f2(G + q,G' + q;r) 
= 2iJ2 [wiG - G") < piG - G") > +n,Ve,iG - G")e-(°-°")'/4] 

G" 

X sin [(G + q) X (G" + q)/2] xg"G' (q; r) 

-2i^l^(q-FG") </9(G-G") > (38) 

G" 

Xsin[(G + q) X (G" + q)/2]xG"G'(q;r). 

Thus the equation of motion for Xgg' (q) at specific wavevector q in the Brillouin zone can be expressed in terms of 
the following matrix equation: 

+ iS - A{q) - D{q)W{q)] x(q) = D{q), (39) 

where 
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AGG'(q) = 2i [W(G - G') < p(G - G') > +nhVeh(G - G')e-(°-°')'/4] 

Xsin[(G + q) X (G' + q)/2] 
Dgg- (q) = -2i < piG - G') > sin [(G + q) x (G' + q)/2] (40) 
VFGG'(q) = W{(i+G)Sgg' 

The collective modes are given by the eigenvalues of the matrix A + DW . Note that generally this matrix is not 
Hermite. However, we can still diagonalize this matrix numerically and it turns out that all the eigenvalues are real 
(and symmetric around zero) for this matrix. 

D. Phonon Shake-up Theory 

We now incorporate the effect of collective modes into the photoluminescence spectrum calculation. Recall that we 
need to calculate the quantity 

nij{Xi,X2;iuJn) = - < TraxAr)ciiT)c]a^x^ > e'^'^^dr. 
■Jo 

Working in the lowest Landau level, the equation of motion for TZij{Xi, X2',t) is given exactly by Eq.(25). In the 
derivation of the terms for contribution from the electron-electron interaction and the electron-hole interaction, by 
applying the many-body Hamiltonian given by Eq.(3), we encounter the following correlation function: 

C,,{XiX2;X3Xi;T',T) = -g < T,at^^(r')ax,(r')ax3(r)c,(r)c](0)at,jO) >, (41) 
In the momentum space, this correlation function can be written as: 

a-,(q,q';r',r) = 1 ^ e-^^^(^^+^=),5x„x.-,„ (42) 
Xi X2 



- e-^«i(^^+^^)<5x3,x,+giC,j(^i,^2;^3,^4;r',r) 



g 

X3X4 



Then the equation of motion for Rij{G,ui) may be written as 



^Rij{G, t) = < p{G, t) > SijSir) - eoRijiG, r) (43) 



-nh 



^V;4G')e'°'x°/2-«"/X-(G-G', 



-| Yl yM)C.A-^, q + G)e'Gx'i/2-^V2 

-| E ^e/.(q)a-,(-q, G)e-''i 
q 

In Eq.(43), 14 (q) and Veh{^) are the Fourier transforms of the electron-electron and electron-hole interactions, respec- 
tively, fo is the energy of the localized hole, and the sum over G' is only over reciprocal lattice vectors, while the sums 
over q are over all wavevectors. R,- specifies the position of the hole in the i th unit cell, and < p{G, r) > e"*^ /■* is 
the expectation value of a Fourier component of the electron density. While this quantity is independent of r in the 
groundstate, it will be convenient for later purposes to formally leave it as an argument of the density. The method 
for computing these Fourier components has been described previously^'^^. Eq.(43) represents the first in an infinite 
series of equations relating an n particle Green's function to the n -\- 1 particle Green's function^®. In the mean-field 
approximation, it was simplified by employing a Hartree-Fock (HE) decomposition of Eq.(41), which converts Eq.(43) 
into a self-consistent equation for Rij' . 

To include shakeup effects, we instead extend this hierarchy to one more level, writing down a self-consistent form for 
Cij which explicitly contains the collective mode excitations. To carry out this program, it is convenient to implicitly 
define a self-energy by writing the last two terms in Eq.(43) as 

-Y,f dr'^iG, G'; r - t')R,,{G', r'), 
^, Jo 
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which implicitly defines the self-energy for this problem. The HF approximation for the PL is equivalent to replacing 
S(G, G'; r - r') by S^^(G, G')S(t - r'), where 

S^^(G, G') = 1^(G - G') < p{G - G'; r) > e'^^^'/^ (44) 

+ :^J2^ehi(l) < Pi-(l,r) > e-^^'/^-"^^^'SG,G-, (45) 
Ztt — ' 
q 

and W is the sum of the direct and exchange Coulomb potentials^' (Note here we keep the negative sign in front 
of q because it is not necessarily a reciprocal lattice vector.) Thus we need to calculate the self-energy 

S = S^^(5(r - r') + (5S 

beyond mean-field approximation 

lation for Cjj 

potential \J which contributes an extra term in the Hamiltonian 



To generate a self-consistent equation for Cjj, we take a functional derivative^® of Eq.(43) with respect to an external 



In doing this, it must be noted that a term directly coupling the Green's function Rij to the external potential U 

-^i?,,(p,r)[/(G-p,r)e^GxP 
p 

must be added to Eq.(43), and the sums over reciprocal lattice vectors must be extended to all wavevectors, because 
\J is not in general commensurate with the lattice. 

We introduce a generalized Green's function F satisfying 

i?,,(p,r) = ^^^(p,p',r) < p(p') > . (46) 
p' 

The equation of motion for the generalized Green's function can be obtained following Eq.(43) 



— _F(p,Pi,r - Ti) = (5ppi(5(r - Ti) - eo-F(p,Pi,-r- n) (47) 



d_ 

-nn J2 V;/.(q)e''ixP/2-5'/4i?(p _ q, p^, r - n 



q 

-V/ rffS(p,q,r-f)^^(q,pi,f-ri), (48) 
q -^0 

Under Hartree-Fock approximation, the equation of motion can be simplified as 

^F"^iP, Pi, r - n) = SppJir - n) - eoF^^ip, pi, r - n) (49) 

-n, J2 V;4q)e''>'P/'"''/'^^''^(P - q, Pi, r - n) 



Comparing Eq.(49) with the equation of motion for the single particle Green's function (Eqs.(12) and (13)), one can 
derive the following relation: 

i^''^(pi,P2;ri - T2) = -e(^°-°-''^)(— =)+'P^xp=/2G(pi -P2;ri -r2). (50) 

Here pi — p2 has to be one of the reciprocal lattice vectors. Using these definitions, we find Eqs. (43) and (48) may 
be rewritten as 
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i^(Pi,P2;ri-r2) = ^^^^(pi,P2;ri -r2) (51) 

-J2 j j dTidhF"^{Y>i,q,i;Ti-Ti) 
qiq2 

X (5i;(qi, q2; fi - -r2)-F(q2,P2;r2 - T2) 

and 

Rij{G,T) = RfP{G,T)- I I dhdhF''^{G,Gi;T-n) (52) 

X ,5S(Gi,G2;n -f2)i?ii(G2,f2). 

Finally, we define the zeroth order Green's functions R^l- and F''-' which satisfy Eqs. (48) and (46) with the self-energy 
S (and external potential) set to zero. It is not difficult to show that 

i?,,(q, r) = i?o^.(q, r)- ^ / / df,df2F°{<i, qi; r - n) (53) 
qiq2"^ 

X i;(qi,q2;-ri - -r2)i?ij (q2 , r2) . 



qiq2 



X [/(qi - q2;r)e''iix'i=/2i?,^.(q2,f). (54) 

It is this form of the equation of motion for Rij that is most convenient for taking a functional derivative. Using 
standard methods^®, we find the Green's functions Rij and Cij satisfy the relation 

^Ig^ = a-,(p', p; r', r) - <,i?,,(p, r) < p(p') 

= <5aj(p',p;r',r). (55) 

Taking the functional derivative of Eq. (54), and defining Cpj^ as the Hartree-Fock decomposition of Cij that turns 
Eq. (43) into Eq. (27), we arrive after much algebra at the relation 



6Cij — SCfj^ + Cfj^ , 



where 



-Rijiq-jyh)- (56) 



(5S(qi,q2)f,-r, 
SUip',r') '-'^^ 

Cfj^ is the correction to Cij beyond mean-field theory, and it will contain the phonon (i.e., collective mode) contri- 
bution to Cij. 

To this point we have made no approximations. To obtain a closed expression for Cfj^ , we need an approximate 
form for S. The best choice available is the Hartree-Fock form, Eq. (45), and so we make this substitution. Combining 
Eqs. (31) and (44), we get: 

,5Sg^(qi,q2)f 
SU{p',T') 

= -l^(qi - q2)e''i^ x'i=/\(p', q2 - qi ; r' - f) 

- 1^ E ^e.(q)e-^^/^-''i^^-<5q,q,x(p', q; r'-r). (57) 
q 

Upon substituting the resulting Cf^^ into the following equation (the phonon mode contribution to Eq.(43)) 



10 



REVTeX 3.0 



Das Sarma Group Preprint, 1996 



J2 j rfr2<5S(Gi,G2;ri-r2)i?ij(G2,r2) 
^^P^(q)e'°^x'l/2C7,^^(-q,Gi + q;n) 



+ I E ^e/.(q)e-^'/4-''i^^'C,^^(-q, Gi; n), (58) 
q 

one obtains the self-energy contributed by collective modes: 

,5S(Gi,G2,ri-r2) 

j2 iqx(Gi-G2)/2-iGixG'i/2+iG2XG^/2 



47r 



;l^(G'i+q)l^(G'2 + q) 

(^^^^(Gi - G'l - q, G2 - G'2 - q; n - r2)xG;G^(q; n - r2) 



W ^ / '''' 

\'^^) r-'r" "'-BZ 



g-(G;+q) = /4-(G^+q) = /4 



xV;;,(G'i + q)V;;,(G'2 + q) 

x^^^^(Gi,G2;ri -r2)xG;G^(q;ri -r2). (59) 

In Eq.(59), Jg^d'^fi represents an integral over wavevectors in the first Brillouin zone of the superlattice. With 
this expression, we are now able to compute the PL intensity. We substitute S = 'S^^S[t — r') + SI] into Eq.(43), 
and Fourier transform this with respect to imaginary time. This means that a Fourier transform of Eq.(59) will be 
necessary, leading to frequency summations of the form (suppressing wavevector arguments) _F(iw„)x(q; — 

ioj„)- To accomplish this, we represent x as a sum over its collective mode poles"'^^; the frequency sums may be then 
computed using standard methods^®. The computation of STi is clearly the bottleneck in this computation, since it 
requires two reciprocal lattice sums and an approximate sum over the continuous wavevector q. We have accomplished 
this using 469 q points in the first Brillouin zone, for one and three electrons per unit cell, which give very similar 
results. Finally, once we have computed ST,, it is straightforward to substitute this into the frequency version of 
Eq.(43), obtain i?(G,w), and from there compute the PL spectrum. 



III. NUMERICAL RESULTS 



Following the algorithm in the proceeding section, one can calculate the collective modes distribution (density of 
states) and their effects on the photoluminescence spectrum beyond the mean-field limit. In this section, we present 
our numerical results for the cases with different strength of electron-hole interaction and different filling factors. We 
will also discuss the temperature dependence of the phonon shakeup effect in photoluminescence from strong magnetic 
field induced Wigner crystal. 



A. General Result 



Shown in Fig. 2(b) is our calculated photoluminescence spectrum for filling fraction i> = 1/5 at T = O(electron 
density Ng = 6 x I0^°cm~^). Our hole is assumed to be strongly localized, and located 250A from the electron 
plane. We ignore the electron-hole interaction in this case. A well-defined shakeup peak may be seen approximately 2 
meV below the main PL peak; a second very weak satellite is observed approximately 3.5meV below the main peak. 
The origins of these peaks may be understood in terms of the phonon DOS, which is illustrated in Fig. 2(a). A van 
Hove singularity, arising from zone-edge phonons, appears as a strong double peak near 0.4meV. Two other peaks 
may be seen near 1.2 meV and 1.9 meV. There are weak sidebands associated with each of these peaks in the PL 
spectrum. The precise interpretation of these peaks is unclear; however, it has been speculated that these represent 
vacancy-interstitial excitations^^. We point out that it is crucial to use a fully quantum mechanical treatment of the 
collective excitations of the lattice to observe these higher order satellites; classical treatments of the phonons^ '"'^"'^ do 
not produce these unusual excitations. 
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(a) v=1 /5 




3.0 



-15.0 -13.0 -11.0 

FIG. 2. Phonon density of states(a) and shakeup spec- 
tmm(b) for a perfect WC with no electron-hole interaction 
for = 1/5. 

It should also be noted that the splitting between the 
main PL peak and the first sideband is actually larger 
than the energy at which the van Hove singularity in 
the phonon DOS appears. The reason for this is that 
there is a strong self-energy renormalization due to the 
phonons in the main PL peak. Physically, this arises be- 
cause the final state of the crystal contains a vacancy, 
which is lowered in energy by a distortion of the lattice 
- i.e., by allowing the electrons surrounding the vacancy 
to relax inward. The self-energy shift accounts for this 
lowering in energy of the final state of the WC, and leads 
to an upward shift of the main PL peak. Our calculation 
shows that the final states in which phonons are excited 
are not nearly so strongly renormalized by lattice relax- 
ation effects, leading to the increased splitting between 
the main PL peak and the sidebands. 



Das Sarma Group Preprint, 1996 

Since in this case, we include a finite number of electrons 
in one supercell to account the lattice distortion due to 
the attraction between the localized hole and electrons, 
there is more structure in the phonon density states as 
shown in Fig. 3(a). 

(a) l/= //5 

50 I ' ' 1 




3.0 



(b) v=l/5 

400 I ' 



'6 






-16.0 -14.0 -12.0 

FIG. 3. Phonon density of states(a) and shakeup spec- 
trum(b) for a WC in the prresence of electron-hole interaction 
for = 1/5. There are 3 electrons per unit cell. 

We also notice that the oscillator strength of the 
shakeup peak increases significantly from the case with 
no electron-hole interaction. Since zone-edge phonons 
correspond to short-range collective modes near the cen- 
ter of the supercell where the hole is located, strong 
electron-hole interaction does enhance their contribution 
to the shakeup spectrum. Some of these local collec- 
tive modes even persist above the melting temperature 
because of the non-uniform electron density distribution 
near the hole. 



B. Effect of Electron-hole Interaction 

Presented in Fig. 3 are the corresponding phonon den- 
sity of states (a) and shake-up spectrum in the presence 
of strong electron-hole interaction for = 1/5 at T = 0. 



C. Other Filling Factors 

Fig. 4 illustrates our results for = 2/7 in the absence 
of electron-hole interactions. We can see splittings in 
the phonon density of states in Fig. 4(a) which are re- 
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lated to the splittings in the Hofstadter spectrum. How- 
ever, zone-edge phonon modes still dominate the shakeup 
spectrum, and as can be seen in Fig. 4(b), except for 
a change in energy scale caused by changing the mag- 
netic field, the photoluminescence lineshape is essentially 
identical to the case of = 1/5. This contrasts sharply 
with the results found in the mean-field approximation^, 
where without electron-hole interactions, a filling v = p/q 
generally yields p distinct lines for a localized hole (see 
Fig.l). While the splittings are so small in that situation 
that they are difficult in practice to resolve, evidently 
shakeup effects wipe out this structure even in principle. 



5- 



(a) v=2/7 




3.0 



600 



t' 400 

s 



(b) v=2/7 



200 



15.0 -13.0 

E-E^-Ef^(raeV) 



11 .0 



FIG. 4. Phonon density of states(a) and shakeup spec- 
tmm(b) for a perfect WC with no electron-hole interaction 
for 1/ = 2/7. 

We note that for the case of an itinerant hole, we expect 
these characteristic splittings to survive shakeup effects. 
The reason is that (neglecting excitonic effects), the PL 
spectrum is related to the product of Green's functions 
for the hole and the electron. For localized holes, only 
the latter has poles in the form of a Hofstadter spectrum 
at the mean-field level, which are wiped out by shakeup 
effects. However, for itinerant holes, the hole Green's 
function also has poles of the Hofstadter form, which are 
unaffected by shakeup. Thus, one should in principle 



be able to identify this characteristic spectrum of bands 
and gaps that is unique to a WC in a magnetic field in 
itinerant hole experiments. 



D. Temperature Dependence 

Fig. 5 illustrate the evolution of the photoluminescence 
spectrum with increasing temperature in the absence of 
electron-hole interaction. As shown in the figure, in the 
solid phase, the main peak slightly shifts to higher en- 
ergy with increasing temperature as in the mean-field 
case. It undergoes a sudden shift upon melting and stays 
at higher energy in the liquid phase. In contrast, the 
shakeup sidebands slightly shift to lower energy with in- 
creasing temperature in the solid phase widening the gap 
between the main peak and the sidebands. However the 
sidebands also shifts rapidly back to higher energy upon 
melting and disappear in the liquid phase. This is nec- 
essarily so, because in the melted phase, the density is 
uniform without electron-hole interaction, and there are 
no collective modes in the lowest Landau level^^. In this 
situation, phonon satellites uniquely distinguish between 
a liquid and a solid state. 



>3 



4000 



3000 



= 1/5 



2000 



^ 1000 




16.0 -12.0 -8.0 



-4.0 



FIG. 5. Temperature dependence of the PL shakeup spec- 
trum for a perfect WC with no electron-hole interaction for 
v = 1/5. Starting from the lowest curve, the correspond- 
ing temperatures comparing to the melting temperature are 
T/Tmeit =0.004, 0.112, 0.224, 0.336, 0.448, 0.561, 0.673, 
0.785, 0.897, 1.009, 1.121, 1.233, 1.345. 

Fig. 6 illustrate the temperature dependence of the 
photoluminescence spectrum in the presence of strong 
electron-hole interaction. In the solid phase, one can see 
that the qualitative behavior of the main peak and the 
sidebands are the same as in the absence of electron- 
hole interaction. However, upon transition, the phonon 
shakeup sidebands undergo a sudden shift to higher en- 
ergy and still persist in the liquid state. In this case, since 
there is a non-uniform electron density near the hole, al- 
lowing some local collective modes to persist even above 
the melting temperature. Except the sudden shifts of the 
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phonon sidebands upon melting, they could not provide 
an definitive evidence for Wigner crystal. However, with 
further increase in temperature, or increased setback be- 
tween the hole and the 2DEG (weakening the electron- 
hole interaction), the oscillator strength of sidebands in 
the liquid phase will significantly decrease. Observation 
of this behavior would distinquish the liquid from the 
solid phase. 



3000 



= 1/5 




14.0 -10.0 -6.0 

E-E^-EjmeV) 



2.0 



FIG. 6. Temperature dependence of the PL shakeup spec- 
trum for a WC in the presence of electron-hole interaction 
for v = 1/5. Starting from the lowest curve, the corre- 
sponding temperatures comparing to the melting temperature 
are T/Tmeit =0.004, 0.045, 0.090, 0.134, 0.224, 0.448, 0.673, 
0.897, 1.121. 

Our results near the melting temperature should be 
viewed with caution, and taken as a qualitative rather 
than a quantitative picture. This is because the melt- 
ing temperature is still calculated within the mean-field 
Hartree-Fock approximation, in which only the contri- 
bution from particle-hole pair excitations is included, 
and thermally excited collective modes are neglected. 
The calculated melting temperature is thus significantly 
higher than the experimentally observed transition tem- 
perature (about 1.5 A'). However, our results should be 
more reliable at low temperatures since the collective 
mode contribution is expected to be insignificant there. 



IV. CONCLUSIONS 

In summary, we have developed a method by which 
shakeup effects in the PL spectrum of a WC from lo- 
calized holes may be computed, that treats the tunnel- 
ing electron and the lattice electrons on an equal foot- 
ing, and uses a fully quantum treatment of the collec- 
tive modes that is realistic over the entire Brillouin zone. 
Our method is quite general, and should be applicable 
to other shakeup problems where quantum fiuctuations 
are important. We find that the Hofstadter spectrum 
found in a mean-field analysis of this experiment is lost 
(although we expect it to survive in itinerant hole exper- 



iments) , and is replaced by a series of sidebands due to 
creation of phonons and other collective excitations of the 
WC. These sidebands are a unique signature of the WC 
which are associated with zone-edge phonons, and can 
in principle be used to distinguish between a liquid and 
crystal state of the electrons. We find that there is a sud- 
den shift in the PL spectrum upon melting of the crystal 
and disappearance of the sidebands in the liquid state in 
the case with weak electron-hole interaction. However, in 
the presence of strong electron-hole interaction sidebands 
persist above melting temperature. 
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